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Abstract. Let A c M. d , d > 2, be a compact convex set and let \i = qq dx be 
a probability measure on A equivalent to the restriction of Lebesgue measure. 
Let v = Qi dx be a probability measure on B r := {x: \x\ < r} equivalent to the 
restriction of Lebesgue measure. We prove that there exists a mapping T such 
that v = jioT^ 1 and T = (f-n, where (p: A — > [0, r] is a continuous potential with 
convex sub-level sets and n is the Gauss map of the corresponding level sets of <p. 
Moreover, T is invertible and essentially unique. Our proof employs the optimal 
transportation techniques. We show that in the case of smooth <p the level sets 
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of ip are governed by the Gauss curvature flow x(s) = — s d_1 K(x) ■ n(x), 

where K is the Gauss curvature. As a by-product one can reprove the existence 
of weak solutions to the classical Gauss curvature flow starting from a convex 
hypersurface. 
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1. Introduction 

vJq The goal of this paper is to introduce a new class of transformations of measures on 

. R d which (heuristically) have the form T = <p ■ Vy/| Vy?| with some function (p. Our 
work is motivated by two intensively developing areas: optimal transportation and 
curvature flows, and establishes an interesting link between these areas. Optimal 
transportation can be described as a problem of optimization of a certain functional 
associated with a pair of measures. The quadratic transportation cost W 2 2 (fi, v) 
between two probability measures /i, v on M d is defined as the minimum of the 
following functional (the Kantorovich functional): 

X 

H| m i— ► / \x\ — x 2 \ 2 dm(xi, x 2 ), mGP(//,v), (1.1) 

JR d xR d 

where "P(//, v) is the set of all probability measures on IR d x IR d with the marginals 
fi and u; here \v\ denotes the Euclidean norm of v G M. d . The problem of mini- 
mizing (11. ip is called the mass transportation problem. This formulation is due to 
Kantorovich [13]. A detailed discussion of the mass transportation problem in this 
setting can be found in [TH]. In many cases there exists a mapping T: IR d — > M. d , 
called the optimal transport between \i and v (or a solution to the Monge problem), 
such that u = fj, o T -1 and 



W%((i,v)= / \x -T(x)\ 2 fx{dx). 

The minimization of the latter integral in the class of measurable mappings T such 
that o T _1 = v is called the Monge problem. If T is a solution to the Monge 
problem, then the image of \x under the mapping x \— > (x, T(x)) to M. d x M. d minimizes 



the Kantorovich functional. However, it may happen that the Monge problem has no 
solution, while the Kantorovich problem is always solvable. It is worth mentioning 
that the first rigorous results related to existence of optimal mappings were obtained 
in the classical work of A.D. Alexandroff [1] on convex surfaces with prescribed 
curvature! If /J and v are absolutely continuous, then, as show by Brenier [B] and 
McCann [16], there exists an optimal transportation T which takes /i to v. Moreover, 
this mapping is /i-unique and has the form T = VW, where W is convex. Under 
broad assumptions, W solves the following nonlinear PDE (the Monge-Ampere 
equation): 

g u (VW) detDlW = Qn, 

where g^ and g u are densities of // and v and D^W is the absolutely continuous 
part of the distributional derivative of D 2 W. At present the optimal transportation 
theory attracts attention of researchers from the most diverse fields, including proba- 
bility, partial differential equations, geometry, and infinite-dimensional analysis (see 
Villani's book [23] and papers [2] and [2TJ). 

The study of curvature flows is a very popular subject in geometry. The theory 
of Ricci flows attracted particular interest after the famous works of G. Perelman 
on the Poincare conjecture. The theory of geometrical flows began, however, with 
flows of embedded manifolds. Let F : M d ~ x — > ~R d be a smooth embedding of a 
smooth compact Riemannian manifold M d ~ x (without boundary). Denote by A the 
enclosed body: dA = Fq(M). We say that an evolution 

F(-, ■): M x [0, T) — > R d 

is a geometrical flow if Fq = F(-,0) and 

^F(x, t) = -g(F(x, t)) ■ n(F(x, t)), (1.2) 

where g : M — > R is some curvature function and n is the outer unit normal vector. 
If g = H is the mean curvature, then F is called the mean curvature flow. If g = K 
is the Gauss curvature, then F is called the Gauss curvature flow. The simplest case 
of the curvature flow is given by the following planar flow: 

= -k{x) ■ n(x), x G M 2 , M = S 1 . 

Consider the flow of closed curves t i— > x(t, •). Under this flow the enclosed volume 
decreases with the constant speed — 2tt. In addition, any non-convex curve becomes 
convex in finite time and then remains convex. Finally, any curve shrinks to a point 
in finite time; the shape of any curve becomes more and more rotund (see [10J and 
|12j). Any multi-dimensional mean curvature flow or Gauss curvature flow starting 
from a convex surface preserves convexity and shrinks the surface to a point (see 
[T5] and [23])- In [S] and [7], equation (11.21) in the case of the mean curvature was 
investigated from the PDE's point of view. It turns out that the surfaces driven by 
(11.21) can be obtained as level sets of a function u(t,x) which satisfies a nonlinear 
degenerate second order parabolic equation of the Monge-Ampere type. A solution 
of this equation is in general understood in some weak sense (viscosity solutions). 
For the PDE approach and viscosity solutions, see the recent book [11]. Concerning 
Gauss flows, see [3]. 

In this paper, we establish the existence and uniqueness of a special measure 
transportation mapping between two probability measures \x and v. It has the 
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following heuristic expression: 

T = tp 



The potential <p has convex sub-level sets. Note that this transportation mapping 
may not be a gradient. Nevertheless, T can be obtained as a degenerate limit 
of some transportation mappings which are constructed by means of the optimal 
transportation techniques. The limiting potential satisfies a degenerate Monge- 
Ampere equation (see the proof of the main theorem). In addition, we show that 
the resulting limit is naturally connected with the Gauss flow. The level sets of the 
potential ip can be associated to a special Gauss flow associated to the measures \i 
and v according to 

x(s) = -s d ^ 1 ^i^. K (x)-n(x). 

Qq{X) 

In the case Qx{x) = \ x \/-i and Qq{x) = we obtain a weak solution to 

x(s) = — c K(x) ■ n(x), 

which is the classical Gauss flow starting from some initial convex hypersurface. 

Finally, we note that in [15], [17], [22] the reader can find other interesting links 
between mass transportation and geometrical flows (in particular, the Ricci flows). 
Some analogs of the presented results in the case of a manifold will be considered 
in our forthcoming joint paper with F.-Y. Wang. 

2. Main result 

Throughout we assume that d > 2 and denote by 7i n the n-dimensional Hausdorff 
measure. For Lebesgue measure we often use another common notation dx. Let 
Int A denote the interior of a set A. 

Recall that, given a compact smooth orientable (d — l)-dimensional surface M 
in M. d , one has the Gauss map n: M — > S" , where n(x) is the global unit outer 
normal vector field. Let Dn: TM r — > TS d 7\ be the differential of n. Choose an 
orthonormal basis {e2, ■ ■ ■ ,e<i} C TM X . Then the matrix Dn can be written as 
(d ei n,ej), where d £i n are the usual partial derivatives of n. The determinant of 
Dn(x) is called the Gauss curvature and is denoted throughout by K(x). 

Below we deal with the case where M is a surface (possibly, non-smooth) of the 
form M = dV, where V is a convex compact set. In this case the normal n(x) is 
well-defined almost everywhere on M. More precisely, for an arbitrary point x G M, 
let us set 

N M>X := { V E S d ~ l ; \/z G V, (rj, z - x) < 0}. 

If Nm, x contains a single element n(x), then n(x) is the unit normal in the usual 
sense. We shall use the fact that one has T-i d ~ 1 {S) = 0, where 

S = {x: Nm, x contains more than one element}. 

Hence the Gauss map n(x) is well-defined 7i d_1 -almost everywhere on M. Moreover, 
one can show that K(x) is well-defined 7Y d_1 -almost everywhere on M (but this fact 
is not used below). 

We shall consider the following Hausdorff distance between nonempty compact 
sets: 

dist(_Bx, B 2 ) = max (sup dist(x, B 2 ), sup dist(x, Bi)j . 

x£B 1 x£B 2 
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Theorem 2.1. Let A C M. d be a compact convex set and let p — g^dx be a probability 
measure on A equivalent to the restriction of Lebesgue measure. Let v = Qi dx be a 
probability measure on B r = {x: \x\ < r} equivalent to the restriction of Lebesgue 
measure. Then, there exist a Borel mapping T : A — > B r and a continuous function 
(p: A — > [0,r] with convex sub-level sets A s = {p < s} such that v — p, o T -1 and 

T — (p • n 7i d -almost everywhere, 

where n = n(x) is a unit outer normal vector to the level set {y: p(y) = p{x)} at 
the point x. 

If p is smooth, the level sets of (p are moving according to the following Gauss 
curvature flow equation: 

x(s) = -3 d - 1 ^^-K{x) ■ n(x) (2.1) 

Q0{X) 

where x(s) G dA r _ s , < s < r, x(0) G OA is any initial point. 

To prove this theorem we develop an approach based on the optimal transporta- 
tion techniques. For every t > 0, we consider a mapping T t that takes p, to v and 
maximizes the functional 



{x,F{x))\F(x)\ t p{dx) (2.2) 
in the class of mappings F with poF^ 1 = v. Equivalently, it minimizes the functional 

Fh/ \x - F{x)\F{x)\ t \ 2 n{dx) 



in the class of mappings F with p o F~ l = v. For t = (12.21) becomes the classical 
Monge-Kantorovich problem. For t ^ standard arguments from the Monge- 
Kantorovich theory show that the set 



{(^(aOITtfc)!*), xe A 



is cyclically monotone, hence belongs to the graph of the gradient of some convex 
function Wt (see [2^ Chapter 2]). This can be shown, for instance, by a cyclical per- 
mutation of small balls (see [21]). More formally, this can be obtained by variation 
of the corresponding Lagrange functional (see [8]). 

If the reader does not want to be concerned with the cyclical monotonicity or 
calculus of variations, we note that VW t is just the optimal transportation of p to 
v o S^ 1 , where S t (x) = x\x\ f . This can be taken for a definition of W t . 

One has the following relations: 



|w t |— 



-, VW t (x)=T t (x)\T t (x)\ t . 



Clearly, |T 4 (x)| < r since T t transforms p into v. 

Throughout the paper we choose Wt in such a way that min^gyi Wt(x) = 0. Define 
a new potential function tp t by 

w t = —p f t +2 . 

t + 2 yt 

One has 



We show below that the limits 

lim cp t = tp, lim T t — T 

exist almost everywhere (for a suitable sequence t n — > oo) and then we prove that 
T is the desired mapping. 

Lemma 2.2. One has 

<Pt < (2 + t)5T* (diam(A))^r5rl, 
f \V<p t {x)\dx< / tp t dH^ 1 < (2 + t)*& (diam(A)) ^r^i H^idA). 

J A JdA 

Proof. By the convexity of W t we have 

W t (x) -W t (y) < (x -y,VW t (x)). 

Choosing y in such a way that W t (yo) = 0, we find W t (x) < diam(A)\VW t (x)\ for 
every x £ A. Since W t = 2^y?t +2 , we obtain 

ft < (t + 2)diam(A)|V^|. 

Let a = . Note that = hence (ptlVptl^ = < r - Therefore, one has 

<Pt = <p]- a <P?< ((2 + t)diam(A)) 1 - a |V^| 1 -V t " 
< ((2 + t)diam(A)) 1 ~ a L|V^ t 



l-a 



< ((2 + t)diam(A)) 1 V = ((2 + t)diam(A)) 2 + l r^. 



By using the convexity of W t again, we get 



o< div( I!^Wdiv^' 



VW t \J V|V^ 

where under div ^ | vv^'l ^) we understand the distributional derivative of the vector 
field j^jHj . Integrating with respect to ip t dx we obtain 

0< / div(^-) Vt dx = - [ \V<p t \dx + I Vt (n A ,^-)dH d - 1 . 
J A MVW J A JdA \ I Vy?t| / 



Hence 



/ \Vip t \dx< [ iptdH' 1 - 1 . 

J A JdA 



Applying the above uniform estimate for (p t we complete the proof. In fact, one 
could do these calculations in the case of smooth densities, where tft has a better 
regularity, and then approximate our densities by smooth ones (the corresponding 
optimal transports converge to VW t ). □ 

Corollary 2.3. There exists a sequence {t n } — > oo such that {tpt n } converges almost 
everywhere to a finite function (p. 

Proof. By Lemma 12.21 every sequence {<pt n } is bounded in W 1 ' 1 (A). By the com- 
pactness of the embedding W l,1 (B) C L l (B) for any ball B C A and the diagonal 
argument (in fact, since in the present situation A is convex, the embedding of the 
whole space H /1 ' 1 (A) is compact), we obtain the claim. □ 
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Lemma 2.4. There exists a sequence {t n } — > oo such that 

lim iV^nl^™ — 1 
almost everywhere, where tp n := (p tn . 

Proof. By Lemma [2.21 one has (p t < (t + 2)diam(A)|V<^|. Hence 

C t \T t \£i < \V<p t \*&, 

where Cf 1 = (2 + t) ~ l (diam(A)) 2+1 . Changing variables one gets the following 
estimate for any 5 > 0: 



Hence 



//(l-C^lV^I^ >5) -0. 



This implies that (1 — C tn \ V<p n \ 1+tn ) tends to zero in fi- measure as t n — > oo. Pass- 
ing to an almost everywhere convergent subsequence one can assume additionally 
that 

lim^jV^I^ > 1 (2.3) 

almost everywhere. Since sup t || V<ft\\L l (dx) < oo by Lemma |2~2| we see that the 

sequence {| V</? n | 1+t « } is bounded in L P (A) for any p < oo. Moreover, by Holder's 
inequality 

Hindoo / \V (p n \ T +^ dx < H d (A). 



A 



Hence, choosing an L p (A)-weak\y convergent subsequence | V</? nm | 1+t «m — > f 7 one 
has 

/ fdx<H d {A). 

J A 

On the other hand, (12. 3ft and Fatou's lemma show that / > 1 a.e., which yields 

lim / |Vy2 n J 1+t ™™ dx—l. 

i 

Hence |Vy? nm | 1+ ' nm — > 1 in the norm of any L P (A), p < oo. Extracting again an 
almost everywhere convergent subsequence we get the claim. □ 
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In what follows we set cp n := (p tn and assume that cp n — > <p and |Vy? n | 1+t ™ — ► 1 
almost everywhere. 

Lemma 2.5. Let C n C B r be convex sets such that Ic n Ic almost everywhere. If 
C is of positive measure, then dist(<9C n , dC) — > 0. 

Proof. The set C can be taken convex by letting Ic := lim „Jr7„ ■ We may assume 
that r = 1. It is known and readily verified by induction that every convex set 
U C B r with H d {U) > 5 contains a ball of volume at least Ki(d)5, where K\(d) 
depends only on d. Let B be a ball of radius e > centered at some point x e dU. 
Then 

H d (Ur\B) > K 2 (d)e d 5, 
where K2(d) depends only on d. It follows that, whenever H d (C n ) > T-C d (C)/2, one 
has || I c - I Cn \\v = n d (CAC n ) > 2- 1 K 2 (d)H d (C)dist(dC n ,dC) d . □ 
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Note that due to convexity one has dist(<9C n , dC) = dist(C n , C). 

Lemma 2.6. The sequence of potentials (pt n converges to (p uniformly on A. In 
particular, p> is continuous and has convex sub-level sets A s = {y: <p(y) < s}. 

Proof. Clearly, it is sufficient to prove the claim for a subsequence. As noted above, 

i 

one can assume, in addition, that | V<£> n | 1+tn — > 1 almost everywhere. Let us redefine 
ip as follows: <p := lim „y„,. Then the sub-level sets A s of (p are convex since the 
corresponding sub-level sets A StTl = {x: (f n (x) < s} of <p n are convex. Since |T n | = 

f>n\ V<£> n | 1+tn , we have shown that \T n \ — > ip almost everywhere. Then the equality 
v = fi o T^ 1 yields that the image of \i under the mapping x i— > (p(x) G M + , denoted 
by fi v G ^(M" 1 "), coincides with v\ x \ G "P(IR + ), where v\ x \ is the image of v under the 
mapping x \— > \x\. Due to our assumptions on u, this implies that fi v has a strictly 
increasing continuous distribution function, i.e., 

1) fJ>(A Sl ) < n(A S2 ) whenever s\ < s 2 , 

2) [i{{(p = t}) = for all t G [0, r}. 

Note that 2) implies that J J 4 s n — > Ia s almost everywhere for each s > 0. By 
Lemma [2.51 we have 

dist(&4 Sjn ,&4 s ) -> 0, s > 0. 

Now, given £ > 0, we divide [0, r] by points s±, . . . , sjv with — Sj| < e and take 
S = maxj<jv dist(dA Si , dA Si+1 ). There exists M such that 

&&{dA Sun ,dA Si ) <5/2 

for every i = 1, . . . , N and every n > M. This implies that sup xgA |<^ n (x) — </ , (a ; )| < 
2e for all n > M. Hence p> n — > 9? uniformly. Since <£> n are continuous as powers of 
convex functions, </? is continuous as well. The proof is complete. □ 

Lemma 2.7. Let N x := N dAip{xhX , where A v ( x ) = {y: <p(y) < tp{x)} and 

S := {x £ A: N x contains more than one element}, 

i.e. S is the set of all the points x such that the boundary of the sub-level set 
containing x is not differentiate at x. Then 7i d (S) = 0. 

Proof. First we consider the case d — 2. Fix an orthonormal basis {ei, and 
identify every unit vector n with a G [0, 2n), where a is the angle between e± and n. 
We write 

n := n a := cos a ■ e\ + sin a ■ e 2 . 
The set S is a countable union of the sets 

S m :={x: [p-q,p + q] C N x }, p, q G Q D [0, 2vr). 

If S has a positive measure, then Ti d {S p ^) > for some p,q. If x G S p ^, then we 
have AtpM C {z: (z — x, n p ) < 0}. Note that the line 

l x , P {z) = {z: (z - x,n p ) = 0} 

intersects S Ptq exactly at x. Indeed, otherwise we get two points x, y such that the 
sub-level sets A^ x ) and A^ y ) both intersect l XiP at two different points and belong to 
the same half-plane P with dP = l XyP . Hence neither A^ x ) C A^ y ) nor A^^ C A^^ 
hold, which is impossible. Thus we obtain that l X;P fl S PtQ = {x}. Finally, applying 
Fubini's theorem and disintegrating Lebesgue measure along the lines parallel to 
Zo, P , we obtain Ti. 2 (S Pjq ) = 0, which is a contradiction. So the lemma is proved 
for d = 2. The multi-dimensional case follows by induction and Fubini's theorem. 
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Indeed, fix an orthonormal basis {ei, . . . , e^}. Note that all sections of a convex 
body are convex. Disintegrating S along e, and applying the result for d — 1, we 
obtain that 

Si = {projection of N x on = has more than one element} 

has measure zero. Since S = UiLi the proof is complete. □ 

Proof of Theorem I2.lt According to Lemma 12.61 and Lemma 12.41 we have 

i 

(p n — > (p uniformly and |Vy5 n | 1+t ™ — > 1 almost everywhere. It remains to prove that 

Vv?„/|Vy?„| -> n 

almost everywhere. Let us fix x G A. Since <p n (x) — > <£>(x), one has lA Vn(x) n Ia < x s 
almost everywhere, where 

According to Lemma 12.71 n(x) is well-defined for almost all x. The same holds 
for every V(p n (x)/\Vtp n (x)\. So, without loss of generality we can fix x in the in- 
terior of A such that n(x) and Vip n (x)/\V<p n (x)\ are well-defined. If the vectors 
V(p n (x)/\V(p n (x)\ do not converge to n(x), then, extracting a convergent subse- 
quence from a sequence of unit vectors {Vip n (x)/\ V(p n (x)\}, we obtain a unit vector 
77 7^ n(x). By using convergence lA Vn(x) „ — * Ia^i one can show that (r),z — x) < 
for all zei, i.e., r] £ N x , which contradicts the choice of x. 

It remains to verify the evolution equation for a smooth potential (p. Indeed, let 
us choose an orthonormal basis {e^} at x such that e% = n and every vector e$, 
2 < i < d, belongs to the tangent space of dA t at x. Let us write the change of 
variables formula for T = tp ■ n. Differentiating along n we find 

d n T = d B (p ■ n + (p ■ d Q n. 

Differentiating the identity (n, n) = 1, we see that <9 n n belongs to the tangent space 
of dA t at x. In addition, d n tp = |Vy>|. Next we note that 

d e T = Lp ■ <9 Ei n, (<9 ei n, n) = 0, 1 < i < d. 

Hence 

det DT= |V^ d_1 det((0 ei n,eA). 

Since K = det((<9 6i n, eA, we have det = | V^ly^ -1 -^. Thus one obtains the 
following change of variables formula (the Monge- Ampere equation): 

Q = Q 1 (<p-n)\V<p\<pf- 1 K. 

It remains to note that the level sets dA s are shrinking with the velocity 1/|V</?| in 
the direction of — n. Hence (12 .ip follows from the change of variables formula. The 
proof is complete. 

Example 2.8. Let A be a convex compact set. Set 

where Cd, r — (yj j — p~l) ■ Varying r we can show the existence of a weak solution 

(in the "transportation sense") to the classical Gauss curvature flow which starts 
from dA and satisfies the equation 

x(s) = — c K(x) ■ n(x), 
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where c can be chosen arbitrarily. 

Certainly, a rigorous justification of this formula requires some additional work, 
since we have not proved that cp is differentiable. 

3. Injectivity and uniqueness 

In this section, we prove that T is invertible and essentially unique. 
Recall that the Legendre transform of a convex function W on a convex set A is 
defined by 

W*(y)=sup((x,y)-W(x)). 

Let dW(x) denote the subdifferential of W at x. Recall also the following known 
fact from the theory of convex functions (see, e.g., {201 Theorem 23.5]). 

Lemma 3.1. Let v G dW(x) for every x G [xi,^], where x\ ^ £2 and [^1,^2] = 
{txi + (1 - t)x 2 , t G [0,1]}. Then [21,0:2] C dW*(v). In particular, W* is not 
differentiable at v. 

In addition to the singular set S C A of all points x such that Nqa uX , where 
t = <p(x), contains more than one element, we introduce another set of degeneracy 
of n defined by 

U — {x G A\S: there is x' G dA t , t = <p(x), such that x' 7^ x and n(x) G Nqa uX '}. 

Proposition 3.2. (i) Consider the set C = dA t for some fixed t. Then the set 
n(U fl C) in S^ 1 has H^ 1 -measure zero. 
(ii) The sets T(U) and 

f(S) := |J <p{x) ■ N dA ^ x 

have v-measure zero. 

Proof, (i) It is sufficient to prove our claim locally on C in a small neighborhood O 
of a point xo where n(xo) is unique. We may assume that n(x ) = — the surface 
C fl O is the graph of a convex function W: B C — > M, where B is an open 
ball containing 0, and that W attains minimum at 0. In addition, we may assume 
that dW(B) is a bounded set. We parameterize C fl O in the following way: 

B 3 (xi, . . . , x d -i) h-> {xx, . . . , x d -i, W(x)). 

Since W is Lipschitzian on B, the surface measure 7i d_1 on C fl O corresponds to 
the measure (1 + |VT4^| 2 ) 2 1 on M.^ 1 . The Gauss map n is given by 

n = ; 1 {-d xl W } -d Xd .W, 1). 

This holds for every (xi, . . . , Xd-i) G B such that (xi, . . . , Xd-i, W(x)) £ S fl C fl O. 
The projection of S fl C fl O on coincides with the points of non-differentiability 
of W. 

It is convenient to identify the half-sphere S d ~ l fl {yd < 0} with its projection 
L Ta!_1 on ]R d_1 and n with the mapping n: , vw taking values in II d_1 . Note 

that the surface measure on has the form mw_i := 1 7i d-1 in the local 

chart (on the set where y\ + • • • + < 1) 

(yi, . . . , y d -i) ' ^ (yi, . . . , \j l -V\ 2/d-i ) • 
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Hence we have to show that 



m d ^(Fo\7W(U')) = 0, (3.1) 
where U' is the corresponding projection of U D C and 



F(x) -- 



X 
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The mapping F is smooth and nondegenerate everywhere. Hence in order to prove 
(jOI it suffices to show that U^iyWiJJ')) = 0. Let us set W := +oo outside of B. 
The Legendre transform W* is finite everywhere. By Lemma [3.11 the set VW(U') 
is contained in the set of nondifferentiability of W*, hence has H ci_1 -measure zero. 

(ii) By Fubini's theorem, it suffices to show that for each t the intersection of the 
set T(U) with the sphere of radius t has zero 7i d_1 -measure. By construction, these 
intersection coincide with the sets T(dA t r\U) defined similarly. Therefore, the claim 
for T{U) follows by assertion (i). 

In order to see that the set T(S) has ^-measure zero, we observe that its inter- 
section with the set T(A\S) of full //-measure belongs to T(U), which is clear from 
the definition of U. □ 

Now we can show that T is invertible. 

Corollary 3.3. The mapping T is injective on a set of full fi-measure. Hence there 
exists a measurable mapping T _1 : B r — ► A such that T(T^ 1 (y)) = y for v-almost 
all y and T~ 1 (T(x)) = x for fi-almost all x. 

Proof Since the equality T(xi) = T[x2) may only happen if <p(xi) = ip(x2), i.e., x± 
and X2 belong to the same level set dA t , it follows from our previous considerations 
that T is injective outside the set T _1 (T(S') U T{U)). This set has /z-measure zero 
because the set T(S) U T(U) has //-measure zero by the above proposition. □ 

Theorem 3.4. The mapping T constructed above is unique in the following sense: 
if a measurable mapping T : A^> B r is such that v = fioT^ 1 and T = (^o' n O; where 
(fQ-. A — > [0, r] is a continuous function with convex sub-level sets A t ^ := {y? < t} 
and n is the corresponding Gauss map, then T = T \i-a.e. 

Proof. Let us show that <po(x) < <^(x) for all x G A. This will yield the equality 
ip Q = ip because otherwise there is t such that /^({<^o < 0) > ^({v ^ 0)> which is 
impossible since both sides equal u(B t ). Set 

C t := [x e A: 16 dA tfi n A t }, D t := |x G A: x E dA t \Int(A ti0 )\ , 

U T :=\JC t , V T :=\jD t . 

t>T t>T 

We observe that for every x G dA t \lnt(A tt o) there exists t' > t depending on x such 
that x G C t >. Indeed, if x G dA tj0 fl dA t , then x G C t . Otherwise one has x G <9A',o 
for some t' = t'(x) > t. Since x G dA t , we have x G lnt(A t /). Hence 

V T C U T . (3.2) 

For every Borel set C C A, set 

f{C) := |J <p(x) ■ N, 



xEC 
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Let us show that 

T (C t ) C T(D t ). (3.3) 

Suppose that xq E dA tt o R A t , t = <po(xo), v = n XQj o E NgA tfi ,x - We show that 
v E T(dA t \Int(A ti0 )). Let us consider the support hyperplane L xo>v 1 u to A tfi 
at Xq. If Xq E dA t and v E NgA t! x , the claim is obvious. Otherwise L XQjV splits 
A t in two convex parts A' t and A t . Since L xo>v is a support hyperplane to A t;0 , one 
of these parts, say, and A 4) o are separated by L XOi „. There exists a hyperplane 
L parallel to L XOjV that is supporting to A t and passes through a point si G cL4 t . 
Then v E NdAtw This proves (13.31) . Hence we have 

T (U t )cT(V t ), 0<t<r. (3.4) 

Suppose now that there exists xq such that ^0(^0) > <p( x o)- Then, by the continuity 
of <p and (po, there is r > for which the inclusion in (13.21) is strict and there is 
a neighborhood in U T not intersecting V T . Therefore, fi(U T ) > fi(V T ). Taking into 
account that T is injective on a full measure set, we obtain 

v(T (U T )) = v(To\T (U T ))) > fi(U T ) > fjL(V T ) = u(T(V T )), 

which contradicts (13. 4p because u(T(V T )) = u(T(V T )) according to Corollary 13. 21 □ 

4. Duality 

Now we consider certain duality properties of the potential if. The duality prin- 
ciple of Kantorovich is a powerful tool for investigating the Monge-Kantorovich 
problem. In our case we also have a kind of the duality formula which relates the 
potential if to some function ijj that can be considered as the support function of the 
family of level sets A t . Note that some interesting duality results for the solution 
of the Monge-Kantorovich problem on a sphere with applications to the prescribed 
Gauss curvature problem have been obtained in [IS] . 

For every y E B r we set 

tl>(y)= SU P ( x iV)- (4- 1 ) 

x- <p(x)<\y\ 

Note that the restriction of ip to dB\ y \ coincides with the support function SY of 
A\ y \ = {x: (p(x) < \y\}, where the support function is defined by 

Sa m {v) := sup (v,x). 

x£A M 

Lemma 4.1. For v-almost all y one has 

^(y) = (T- 1 (y),y). (4.2) 

Proof. It is clear that the supremum on the right-hand side of (14.11) is attained at 
a point p such that y E Nqa vM>p - This implies that p coincides with T~ l (y) for v- 
almost all y, hence is 1/- almost everywhere well-defined, which yields our claim. □ 

Now we show how to describe ip as a limit of certain functions depending on pre- 
limit potentials ip t - Recall that the Legendre transform W* satisfies the inequality 

W t (x)+W t *(y)>(x,y). (4.3) 

An equality holds if and only if y E dW t (x) and x E dW*(y). Moreover, W t and 
W£ satisfy the identities 

W; o VW t (x) = x, VW t o VW t *(y) = y 
11 



almost everywhere on the sets A and VWt(A). Since 

VW t = \T t \%, 

one has 

T t -\y) = VW7(|y|*y). 
In what follows we denote by I the identity matrix and by I z the orthogonal projector 
on the one- dimensional vector subspace generated by z, i.e., 

(v,z) z 



Lv 



\z\ \z\ 



We have found a sequence t n — > +00 for which the mappings T tn converge to 
T almost everywhere on A, hence converges in measure fi. For this sequence, the 
following holds. 

Lemma 4.2. The mappings T^ 1 converge to T _1 in measure v . Hence there exists 
a subsequence t' n — > 00 such that T^ 1 — > T _1 v-almost everywhere. 

Proof. Since /x o T -1 =/jo T^ 1 = 2/, for any //-measurable function /, the functions 
/ o T tn converge to / o T in measure /i (see, e.g., [U Corollary 9.9.11]). Hence the 
mappings T~ x oT tn converge to T -1 oT = J in measure //. Therefore, for every c > 
one has 

u(y: \T t -\y)-T-\y)\>c)=^{xeA: \x - T^T^x)] > c) - 
as n — > 00, which completes the proof. □ 
Theorem 4.3. Let a function ijj t be defined by the relation 

Wt(z) = \z\^tjjt(z\z\~^). 

Equivalently, 

My) 



W;{y\y\ 


') 




y\ 


t 



Then one has ip = lim ip tn almost everywhere for some sequence {t n }. 

Proof. Note that it is consistent with our previous choice of Wt to assume that 
W?(0) = 0. Indeed, W t (x) > (x,y) - W£(y), hence taking y = we find W t (x) > 0. 
Taking any x G dW^(0) we easily obtain W t (x ) = 0. Indeed, for (x ,0) inequality 
(14.31) becomes an equality, hence W t (x ) + W t *(0) = (x , 0) = 0. 

The inequality W t *{a) - W t *(b) < (a - b,VW*(a)) yields, by substituting 6 = 
and a = y|y|*, that 

My)< ( yw:(y\y\ t ),y) = (y,T-\y)). 

Similarly, if a = and 6 = y\y\ l , one has ^(y) > (v,y) for any u G cW t *(0). In 
particular, 

|^t(y)| < diam(A)|y|. 

One has 

l + t \z\ V 1 + 1 J 

Substituting z = \y^y we obtain 

T t -\y) = j^Mv)^ + (1- VtMv). (4-4) 
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Taking the inner product with y we find 

(T t -\y),y) = ~^- t My) + j^<Wh{v),v)- ( 4 - 5 ) 

In view of Lemma l4.2l and equality (14. 2p it suffices to show that j^-(Vipt n {y)y y) ~ > 
almost everywhere for some {t n }. Indeed, since ip t < (Tt{y),y), we obtain from 
(USD that 

(Tf\y),y) < Y^(TT\y),y) + j^WMvU)- 

Hence 

il>t<(T t -\v),y)<(VMv),v)- 
Taking into account that ip t > — diam(A)\y\, we see that {(Vip t (y),y)} is uniformly 
bounded from below. The integration by parts formula yields 

{Vil) t {y),y)dx= -d ip t dx + {y^tdH^ 1 . 

Applying again the estimate \ip t \ < diam(v4)||/| we obtain 

sup / (Vip t (y),y) dx < oo, 

t J B r 

hence sup f \\(V^t(y),y)\\L^(B r ) < oo. Therefore, 

Km-^-||(VV t (y),2/>|Ui(B r )=0. 

Extracting a subsequence we complete the proof. □ 

Remark 4.4. Taking a scalar product of (I4.4p with any vector V-Ly we obtain the 
equality d v ip n (y) = (T~ l (y),v). Let us set 

d v tp(y) := lim d v if) n (y). 
In view of convergence T n — > T this definition makes sense. Moreover, we have 

dMy) = (T-\y),v), (4.6) 

for any vl.y. 

Taking into account f!4.2p we obtain the following remarkable relation: 

ip(y) 



T-\y) = TT^) + ^d ei(y) ^(y) ei (y), 



lt71 i=2 

where {ei(y)} is an orthonormal system of unit vectors chosen in such a way that 
ei(y) = y/\y\ and ei(y)±y, 2 < i < d. 

Remark 4.5. Let us see what happens in the limit with the duality formula 

W t (x)+W*(z) > (x,z). 

It can be rewritten as 

' <f t + 2 (x) + \y\ t i/; t (y)>(x,y)\y\ t 



t + 2 
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by letting z := y\y\ l - If y — T t (x) = yj t y^j| V<^| 1+t , then an equality holds. It is 
known that this is possible only if the pair (x, y) belongs to the graph of T t . Hence 
we obtain the following duality relation: 

J-^ t 2 (x)|V^(x)|-^ +1 H (T t (x)) = (x,T t (x)). 

In the limit t — > oo we find 

1>{T{x)) = (x,T(x)). 

It is worth noting that for the constructed transformation T one can establish a 
change of variables formula involving certain analogs of Alexandroff 's determinants. 
Such formulas which neglect singular components are known for optimal transfor- 
mations and triangular transformations (concerning the latter, see [3J Ch. 10] and 

0)- 
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